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Abstract 

We investigate the phase diagram of a three-component system of 
particles on a one-dimensional filled lattice, or equivalently of a one- 
dimensional three-state Potts model, with reflection asymmetric mean 
field interactions. The three types of particles are designated as A, B, 
and C. The system is described by a grand canonical ensemble with 
temperature T and chemical potentials TXa, TXb, and TXc- We find 
that for Xa = Xb = Xc the system undergoes a phase transition from 
a uniform density to a continuum of phases at a critical temperature 
T c = (27r/v / 3)~ 1 - For other values of the chemical potentials the sys- 
tem has a unique equilibrium state. As is the case for the canonical 
ensemble for this ABC model, the grand canonical ensemble is the 
stationary measure satisfying detailed balance for a natural dynam- 
ics. We note that T c = 3T C , where T c is the critical temperature for 
a similar transition in the canonical ensemble at fixed equal densities 
r A = r B = r c = 1/3. 

Keywords: ABC model, grand canonical ensemble, reflection asymmetric 
mean field three-state Potts model 



1 Introduction 



In this paper we study the phase diagram of the three species ABC model 
on an interval as a function of the chemical potentials and the temperature. 
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The system is defined microscopically on a lattice of N sites in which each 
site is occupied by either an A, a B, or a C particle. The energy is of mean 
field type, with an interaction which has cyclic symmetry in the particle 
types but is reflection asymmetric: 

. N-l N 

En (Q = Jj12 lvc(i)VB(j) + VA(i)vc{j) + VB{i)VA{j)}- (1.1) 

i=l .7=1+1 

Here the configuration £ of the model is an iV-tuple (£i, . . . , (n), with Q = A, 
B, or C, and Tj a {i), a = A,B, C, is a random variable which specifies whether 
a particle of species a is present at site i: r] a (i) = 1 if Q = & and r/ a (i) = 
otherwise, so that always Yla^ii ) = 

We remark that we may also regard the model as a reflection asymmetric 
mean field three state Potts model. The asymmetry of the interaction, 
however, gives this model system very different behavior from that of the 
usual symmetric mean field model pp. Similar but short-range (in fact, 
nearest neighbor) reflection asymmetric interactions occur in chiral clock 
models [2J[3l|4]; see Remark [5J] below. 

The equilibrium probability of a configuration £ is given by the grand 
canonical Gibbs measure 

M(S)A (C) = H- 1 exp[-PE N (g + A a JVa(C)], (1-2) 

a 

where f3 is the inverse temperature, Aa, Ab, and Ac are f3 times the chemical 
potentials, N a = YliLiVa(i) with J2 a N a = N, and 5 is the usual grand 
canonical partition function. We prove here that in the scaling limit (N — > 
oo, i/N — > x G [0,1]) the equilibrium density profiles p{x) are unique and 
spatially nonuniform when the A Q 's are not all the same. When Aa = A# = 
Ac the densities are spatially uniform above a critical temperature T c = /3~ , 
with /3 C = 2tt/\^3; below T c the profiles have a natural extension to periodic 
functions with a period three times the length of the system. 

One may compare the behavior described above with that of the same 
system in the canonical ensemble, in which the -/V a are taken as fixed; this is 
the only case considered previously. The results are quite different, that is, 
we have inequivalence of ensembles (see [U [6] for recent reviews) . We give in 
Section 2 a brief history of the ABC model with fixed particle number and 
a summary of results for that system. In Section 3 we describe a stochastic 
evolution satisfying detailed balance with respect to the measures /^s,a(C)) 
and in Section 4 we establish the phase diagram. Section 5 gives a discussion 
of some related models and problems. 
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2 The ABC model in the canonical ensemble 

The ABC model was introduced by Evans et al. [7] (see also [HI l9l flQl fTTl [T2l 
[TBI as a one dimensional system consisting of three species of particles, 
labeled A, B, C, on a ring containing N lattice sites; we will typically let 
a = A, B, or C denote a particle type, and make the convention that 
a + 1, a + 2, . . . denote the particle types which are successors to a in the 
cyclic order ABC. The system evolves by nearest neighbor exchanges with 
asymmetric rates: if sites i and i + 1 are occupied by particles of different 
types a and 7, respectively, then the exchange 07 — > 7 a occurs at rate 
<7<lif7 = a + l and at rate 1 if 7 = a — 1. The total numbers N a of 
particles of each species are conserved and satisfy Yl a N a = N. In the limit 
N — > 00 with N a /N — > r a , where r a > for all a, the system segregates 
into pure A, B, and C regions, with rotationally invariant distribution of 
the phase boundaries. 

In the weakly asymmetric version of the system introduced by Clincy et 
al. [9], in which q = , the stationary state for the equal density case 

Na = Nb = Nc is a Gibbs measure of the form exp{— /3En}, so that the 
parameter ft = T _1 plays the role of an inverse temperature. The energy 
Ej\f is given by (jl.ip , and the condition Na = Nb = Nc ensures that this is 
translation invariant, despite the appearance of a preferred starting site for 
the summations. 

Ayyer et. al. [H] studied the weakly asymmetric system on an interval, 
that is, again on a one-dimensional lattice of N sites but now with zero flux 
boundary conditions, so that a particle at site i = 1 (respectively i = N) can 
only jump to the right (respectively left). For this system the steady state 
is always Gibbsian, given by exp{— {3En} with En as in (jl.ip . whatever the 
values of ./V^, Nb, and Nc- When Na = Nb = Nc the steady state of 
the system thus agrees with that on the ring, so that the invariance under 
rotations on the ring then implies a rather surprising "rotation" invariance 
of the Gibbs state on the interval. We describe the results of [13] in some 
detail, since the work of the current paper depends heavily on them. 

To identify typical coarse-grained density profiles at large N, [2] con- 
siders the scaling limit 

iV^oo, i/N-tx, a; €[0,1]. (2.1) 

For this limit there exists a Helmholtz free energy functional /3~ 1 J 7 ({n}) 
of the density profile n{x) = (nA(x),n,B(x),nc(x)). F is the difference of 
contributions from the energy and entropy: 

T({n})=p£({n})-S({n}), (2.2) 
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where £ {{n}) and S ({n}) are the limiting values of the energy and entropy 
per site: 



£ (W) 

S({n}) 



\ dx dzVn a (x)n a+2 (z), (2.3) 

J Jx a 

— / dx \^ n a (x) In n a (x). (2-4) 
JO 



We will write J- = J-^ when we need to indicate explicitly the /3 depen- 
dence. Only the canonical ensemble was considered in [14J, so that for some 
fixed positive mean densities r^, rs, rc satisfying ta + tb + re = 1 the 
profiles n{x) in (|2.2p - (|2,4p satisfy the conditions 



< n a (x) < 1, 2_, n a{x) = 1, and / 
« Jo 



n a (x)dx = r a . (2.5) 



The typical profiles in the scaling limit are those which minimize T\ it was 
shown in [14J that such minimizers always exist and satisfy the ELE derived 
from T . To obtain the ELE one defines 

T a (x) 



5n a (x) 

logn a (x)+f} I [n a+ i(z) - n a+2 {z)} dz + 1 + (3r a+2 (2.6) 



to be the variational derivative taken as if the profiles ua(x), ub(x), and 
nc(x) were independent; the constraints (j2.5j) then imply that at a station- 
ary point of T both Ta — J~c an d J~b — J~c are constant. After simple 
manipulations (see also Section U] below) this yields the ELE satisfied by the 
typical profiles p(x): 

^T L =Ppa(j>a-l-Pa+l), a = A,B,C. (2.7) 
dx 

These are to be solved subject to (j2.5j) (written in terms of p rather than 
n). 

It follows from (|2.7p that all relevant solutions satisfy naPa( rE ) = K for 
some constant K with < K < 1/27. For K = 1/27 they are constant, 
with value 1/3; for K < 1/27 they have the form 

Pa(x) =y K (2P(x-l/2)+t a ), 0<x<l, (2.8) 
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with yxit) a solution, periodic with period tr-, of the equation 

± y ' 2 + ±Ky-±y 2 (l-y) 2 = 0; (2.9) 

here t = 2f3x + constant, yx is uniquely specified by requiring that it take 
on its minimum value at the points t = titk, n £ Z. The phase shifts t a in 
(l2~8l) satisfy 

tA = t B +T K /3 and t c = t B -r K /3. (2.10) 

Remark 2.1. Equation (j2.9j) describes a particle of unit mass and zero 
energy oscillating in a potential U K {y) = Ky/2-y 2 (l -y) 2 /8. The constant 
solution y = 1/3 appears for IT = 1/27. For K < 1/27, yx(t) is an even 
function which is strictly increasing on the interval [0, rjf/2]; it was shown 
in [TJ] that tk is a strictly decreasing function of K. Because the potential 
is quartic in y the solution is an elliptic function. Further properties of the 
function yx are summarized in Proposition I A. II of Appendix lAl 

Equation (|2.8j) indicates that nonconstant solutions of the ELE are ob- 
tained by viewing yx(t), an d its translates by one-third and two-thirds of a 
period, in some "window" of length 2/3. If one is given (3 and r = (r^, r#, re) 
then one must determine K and one of the phase shifts, say t#, so that 

1 f 13 

r a = 7^ y K (t + t a )dt, a = A,B,C. (2.11) 
Z P J-p 

The solutions which minimize J- were completely determined in [T3]. In 
stating the result, we use the following terminology: a solution is of type n 
if (n — 1)tk < 2/3 < nrx, that is, if the window contains more than n — 1 
and at most n periods of the function yx- 

Theorem 2.2. Suppose that ta, tb, and tq are strictly positive. Then: 

(a) If rA = rs = rc = 1/3 then for the equations (12.T|) with (12. 5p there 
exist (i) the constant solution, (ii) for /3 > nf3 c = 2irn\^3, n = 1,2,..., a 
family of solutions, of period tk = 2/3/re and hence of type n, differing by 
translation, and (Hi) no other solutions. The minimizers of the free energy 
are, for (3 < (3 C , the (unique) constant solution and, for (3 > j3 c , any type 1 
solution. 

(b) For values of r other than (1/3, 1/3, 1/3) there exists for all f3 a unique 
type 1 solution of these equations which is a minimizer of the free energy. 

(c) At zero temperature (/3 — )■ oo) the system segregates into either three or 
four blocks, each containing particles of only one type. 
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3 Dynamics of the grand canonical ABC model 

We now turn to consideration of the ABC model on the interval when the 
number of particles can fluctuate; we will abbreviate this as the GCABC 
model. In Section 1 the corresponding grand canonical measure ug \ (see 
(jl.2p ) was presented in the equilibrium setting as a Gibbs measure obtained 
from the energy function (jl.ip and chemical potentials /3 _1 A a . Just as for 
the canonical Gibbs measure, however, one may alternatively view this as 
the stationary measure for some dynamics; we describe two possibilities here. 
(A different generalization of the ABC model to a nonconserving dynamics, 
in which the system is on a ring, vacancies are permitted, and the total 
number of particles fluctuates but the differences N a — iV 7 are conserved, is 
given in |151 116| . When all the N a are equal the stationary measure has the 
form of a grand canonical ensemble.) 

In the first dynamics we consider there are particle exchanges between 
adjacent sites, with the same rates as for the canonical dynamics. To allow 
the number of particles to fluctuate, however, we introduce two new possible 
transitions. First, if the particle at site i = 1 is of type a then with a rate 
equal to Ce~ Xa the entire configuration is translated by one site to the left, 
the particle at site i = 1 disappears, and a particle of species a + 1 is created 
at site i = N . Second, with a rate equal to Ce~ Xa+1 the reverse transition 
occurs. Here C is a constant which we shall in the future take equal to 1. 
This dynamics satisfies the detailed balance condition with respect to the 
Gibbs measure (jl.2p : if a transition £ — > £ arises from an exchange of 
particles the argument is as for the canonical model [14], while if it comes 
from a transition of the new type, say in the "forward" direction as described 
above, then E^(Q = En(C') but N a decreases by 1 and N a+ i increases by 
1, and the detailed balance condition e~^ a vp{C l ) = e~ Xa+1 u^(^') follows. 

Remark 3.1. One may also obtain this dynamics by considering a ring of N 
sites, with each site occupied by an A, B, or C particle and with a marker 
located on one of the bonds between adjacent sites. Adjacent particles 
exchange across any unmarked bond with the usual ABC rates, while the 
marker may move one bond to its left or right, and in doing so it changes 
the species of the particle it passes: with -x- and — denoting a marked and 
unmarked bond, respectively, the transition -x- a - — > — (a + 1) -x- 
occurs with a rate equal to e~ Xa and the reverse transition with a rate 
equal to e~ Xa+1 . If one then obtains a configuration on the interval from 
a ring configuration by letting the marked bond identify the boundaries of 
the interval — effectively by cutting the ring at the marked bond — one sees 
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easily that the inherited dynamics on the interval is precisely the dynamics 
discussed above. A slight variation of this idea was mentioned in [33] • 

We define the second dynamics only for the case in which all the A Q are 
equal. We obtain it by first defining a dynamics for the constrained ring: a 
ring of 3iV sites populated by A, B, and C particles but with a restriction 
to configurations (&)f=i which satisfy 

i i+N = & + 1 (3.i) 

(addition on the site index is modulo 3JV); that is, if an A particle is on 
site i then there must be a B particle at site i + N and a C particle at site 
i + 2N, etc. The dynamics for the constrained ring is given by a modification 
of the usual rules of the canonical ABC model on a ring: exchanges occur 
simultaneously across three equally spaced, unmarked bonds in the usual 
ABC manner, with rate 1 for the favored exchanges and rate q = e~^l N for 
the unfavored ones. 

We consider now any fixed block of N consecutive sites on the con- 
strained ring and ask for the induced dynamics on configurations in this 
block. Two types of transitions occur: nearest-neighbor exchanges at stan- 
dard ABC rates for a system of size ./V and inverse temperature (3 (i.e., rates 
1 and q = e~@/ N ) and a transition corresponding to an exchange on the con- 
strained ring across the boundaries of the block. To understand the latter, 
suppose the configuration within the block has the form (a + 2) £ (a + 1), 
with £ any configuration on JV-2 sites; then (|3.ip implies that the particles 
immediately to the left and right of the block are of type a, and a transition 
from (a + 2) £ (a + 1) to aQa occurs at rate 1. The reverse transition oc- 
curs at rate q, and no such transition occurs when the block configuration is 
(a + 2) £ a. Then using A a = A# = Ac one checks, just as for the dynamics 
considered above, that if one identifies the block with an interval of N sites 
then this dynamics satisfies the detailed balance condition with respect to 
the grand canonical Gibbs measure (II. 2p . 

On the constrained ring there are equal numbers of A, B, and C particles, 
from (I3.ip . so that the energy E^n (that is, the energy given by fll.lf) with N 
replaced by 3A throughout), and thus the restriction of the Gibbs measure 
Z~ x exp{— (5Ezn} to particle configurations satisfying (13. li . is well defined 
and independent of the starting point of the summations [TJ. Moreover, this 
is the invariant measure for the constrained ring dynamics defined above, as 
one again checks by verifying detailed balance. With the discussion above 
this shows that the restriction of Z -1 exp{— /3E^]y} to the block of N sites 
is the Gibbs measure (ll.2p . One may also verify this from the fact that if £ 
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is a constrained ring configuration and £ the portion of that configuration 
within the block then 

E 3N (0 = E N (Q + N/3. (3.2) 

Thus we can study the GCABC with \a = Xb = Ac* by studying directly 
the constrained ring. 

3.1 The scaling limit for the constrained ring 

To identify typical coarse-grained density profiles at large N on the con- 
strained ring we consider the scaling limit (|2.ip with N replaced by 3N 
(N — > oo with i/3N — > x G [0,1]) and find the appropriate free energy 
functional. The scaling limit of the energy per site is still given by (|2.3p . 
but because the full microscopic configuration under the constraint (]3,ip is 
determined by the configuration of the first N sites the entropy per unit site 
is only 1/3 of (I2.4h . This leads to a free energy functional 

f}£{{n})- l -S({n})= l ~TW){{n}). (3.3) 

Here £({n}) and S({n}) are as in ()2.3[) and (|2.4p and n is a constrained 
density profile, that is, one which satisfies the continuum equivalent of (13. ip : 

n a (x) = n a+ x(x + 1/3), (3.4) 

where the addition x + 1/3 is taken modulo 1. J 7 ^) is the free energy 
functional at temperature 3/3 of the (unconstrained) canonical system on an 
interval, as defined in (|2.4p ; equivalently, because there are equal numbers 
of particles of each species, this is the free energy functional on a ring |14j . 

Typical (coarse-grained) profiles at inverse temperature f3 on the con- 
strained ring, for large N, correspond then to continuum density profiles 
p(x) which satisfy the constraint (13. 4p and minimize the free energy over all 
such constrained profiles. It follows from (|2.8p and ()2.10p . however, that the 
typical profiles (minimizers) for the canonical free energy, which are a priori 
unconstrained, do in fact satisfy (|3.4p . Thus by (|3.3p the typical profiles for 
the constrained ring are the same as the typical profiles of an unconstrained 
system on the ring at inverse temperature 3/3. This is illustrated in Figure [IJ 
where we plot time-averaged profiles from Monte-Carlo simulations of the 
constrained ring at /3 = 10.152 and the exact solution |14j for the uncon- 
strained ring at /3 = 30.456, showing close agreement. (We can use time 
averaging rather than spatial coarse graining for this comparison because 
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Figure 1: Typical profiles in a large system. The dotted curves are time- 
averaged occupation numbers in a constrained ring of size N = 1800 at 
inverse temperature j3 = 10.152. The solid curves are the corresponding 
elliptic functions obtained from the exact solution of [14] at temperature 
P = 30.456. 

the time scale for the profile to drift around the ring is much larger than the 
simulation time scale.) 

It follows from this discussion that when the chemical potentials are 
equal the critical temperature f3 c for the grand canonical ensemble on an 
interval, which is represented by the part of the constrained ring between two 
markers, is f3 c = (3 c /3. Typical configurations are constant if /3 < j3 c and for 
/3 > /3 C are a portion of the typical profile for the canonical system at inverse 
temperature 3/3; the latter is periodic and in the GCABC system we see a 
randomly-selected one-third of a period. These properties are confirmed in 
Section |4] by direct analysis of the grand canonical system in the scaling 
limit. 

4 The phase diagram of the GCABC model 

In this section we discuss the GCABC model directly in the scaling limit 
(|2.ip . From (jl.2p we see that the new free energy functional F({n}) (= 
J 7 p,x({ n })) i which is the negative of the pressure multiplied by j3, is obtained 
by adding chemical potential terms to the free energy functional of the 
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canonical model: 



^({n})=J-({n})-^A Q f 
„ JO 



dx n a (x), 



with T given by ()2.2|) . The profiles now are constrained only by 
< n a (x) < 1 and ^^n Q (x) = 1. 



(4.1) 



(4.2) 



We will always normalize the chemical potentials so that Y2a = (al- 
though with this normalization we cannot conveniently consider the limit in 
which just one of the A Q becomes infinite). Just as for the canonical model 
|14j it can be shown on general grounds that for every j3, A the free energy 
functional has at least one minimizing profile p(x) which belongs to the in- 
terior of the constraint region, i.e., satisfies < p a (x) < 1 for all a,x (and 
of course J2 Q p a (%) = 1 for all x). Prom this it follows that p(x) will satisfy 



5p A (x) 
5 

5p B (x) 



Hip}) 
H{p}) 



PC=1-PA~PB 



PC= l ~PA-PB 



(F A (x) - X A ) - {Tcix) - A c ) = 0, (4.3) 
(T B (x) - X B ) - {Tcix) - X c ) = 0, (4.4) 



with J-" Q (x) as in ()2.6p . so that F a (x) — X a is independent of a. But one 
finds from (|2.6p that p a dT a (x) / dx = 0, so that 



F a (x) - X c 



C 



(4.5) 



for some C independent of x and a. Differentiating (|4.5p leads again to 

^ = /3pa(p a -i ~ Pa+i), a = A,B,C. (4.6) 
ax 



Moreover, (|4.5p implies that J-" Q (0) — X a = 7^+1(1) — A a+ i, which with (|2.6 
yields the boundary condition 



p a (0)e~ 



Pa+l 



(iK 



A,B,C. 



(4.7) 



Note that (|4.7p is consistent with the (first) dynamics described in Sec- 
tion sec:dynamics. 

Equations (j4.6|) and (j4.7|) may be taken as the ELE of the model (it is 
easy to verify that these imply (|4.5p ). Solutions of (|4.6p are, by the analysis 
of |14| , of the form (|2.8p , with phase shifts satisfying (|2.10p . It remains only 
to consider the effect of the boundary condition (|4.7p . 
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Let us begin by considering the case A^ = A# = Ac*, in which (|4.7p 
becomes p a (0) = p a+ \{l). Certainly the constant profile with p a (x) = 1/3 
for all a, x satisfies this condition and hence is a solution for all f3. From 
(|2,10p we see that a nonconstant solution (|2.8p will satisfy this condition if 
and only if 

VK(t a -P) = y K (t a + P-TK/a), a = A,B,C. (4.8) 

The properties of uk mentioned in Remark 12.11 imply that (|4.8p can hold 
if and only if (t a — /3) ± (t a + /3 — r^-/3) is an integer multiple of tk- The 
choice of the positive sign here leads to no solutions consistent with (|2.10p ; 
the negative sign gives 2/3 = (3n — 2)r^/3 for n = 1,2,3, .. .. Since the 
minimal period of solutions of (|2.7p is 2/3 c = 47r\/3, a nonconstant solution 
of (|2.7p and (|4.7p can exist only if /3 > /3 c /3; thus as in Section[3]we find that 
Pc = Pc/3 is the critical inverse temperature for the GCABC model. There 
is no constraint on the t a other than (|2.10p , so that there is a one-parameter 
family of solutions differing by translation. 

Following the usage of |14j it is natural to refer to the solutions just 
discussed for which 2/3 = (3n — 2)r^-/3 as being of type n. We will, again 
as in [2], extend this classification to the case of general A: a solution (|2.8p 
of (I4.6P and (14. 7p will be said to be of type 1 if 2/3 < t/^/3 and of type n, 
n = 2,3,..., if (3n- 5)t k /3 < 2/3 < (3n - 2)t k /3. With this terminology 
we can state our main result. 

Theorem 4.1. (a) If Xa = A# = Ac* then for the equations fj4.6|) and 
(|4.7p there exist (i) the constant solution, (ii) for (3 > (n — 2/3)/3 c = 
2ir(n — 2/3)\/3, n = 1,2,3,..., a family of solutions of type n, differing 
by translation, and (Hi) no other solutions. The minimizers of the free en- 
ergy functional J- are, for ft < ft c /3, the (unique) constant solution and, for 
(3 > p c , any type 1 solution. 

(b) If not all A a are equal then there exists for all P a unique minimizer of 
the free energy functional JF; moreover, this minimizer is a type 1 solution 
of dUED and (glTD - 

We give the proof of part (a) of this theorem here; the more technical 
proof of (b) is presented in Appendix lAl 

Proof of Theorem \4-l^ a): The discussion at the beginning of this section 
establishes the first statement of the theorem; it remains to show that the 
type 1 solution, when it exists, minimizes the free energy. We do so by 
reducing this problem to the corresponding one for the canonical ensemble; 
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the argument is similar to the consideration of the constrained ring system 
in Section [3l For any profile n(x) = (ua{x), tib(x), nc{x)) (where it is 
understood that < n a (x) < 1 and ^ Q n a (x) = 1) define the tripled profile 
0({n}) by 



(e({n})Ux) 



n a (3x), if < x < 1/3, 

n a -i(3x-l), if 1/3 < x < 2/3, (4.9) 
i n Q _ 2 (3x-2), if2/3<x<l. 



The profiles which have the form 0({n}) for some n are precisely those 
satisfying ([3.4P ; in particular, each 0({n}) gives equal mean densities to the 
three species. 

Now a simple computation shows that for any profile {n a (x)}, 

F^\{n}) = ^ 3/3 )({6({n})}) - 0/3. (4.10) 

(Note that this free energy differs by an overall factor, plus an additive 
constant, from that of (I3.3P ; the difference arises because here we started 
from the energy and entropy per site on the interval of size N, and in f|3.3|) 
from the energy and entropy per site on the ring of size 37V.) Thus the 
problem of finding the minimizer(s) of ,F^' )({n}) = J-^({n}) over all 
profiles n is equivalent to finding the minimizer(s) of J-^\{n}) over all 
profiles satisfying (|3.4p . On the other hand, the minimizers of over 
all equal-density profiles are given in Theorem 12. 2f a): the constant solution 
if 3/3 < j3 c and the solution of (minimal) period 6/3 if 3/3 > f3 c (this is the 
type 1 solution for the canonical model). Because these are either constant 
or periodic, they satisfy fj3.4|) and hence are also the minimizers over all 
such profiles. But these minimizers are precisely the images under of the 
profiles identified as minimizers in Theorem 14.1( a). □ 



Remark 4.2. In the argument above the essential role of the tripling map 
is to convert the problem of minimizing T with respect to arbitrary vari- 
ations in the profiles to the previously solved problem of minimizing under 
variations which preserve the condition f Q dxn a (x) = 1/3. Other conclu- 
sions may be obtained similarly; we mention briefly two examples. 

(a) It was shown in [U] that, for /3 < (2/3\/3)/3 c and any r = (ta, ?\b, re), 
J~({n}) is convex as a functional of profiles satisfying (|2,5p . Via this 
implies that J^({n}) is, for (3 < (2/3v / 3)/3 c , convex as a function of profiles 
satisfying (I4.2j) . 
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(b) The two point correlation functions on the interval are related to those 
on the constrained ring by 

(n(x)n(y)) intcrval = (n(x/3)n(y/3)) ring - (4.11) 

The latter (denoted below simply as (•)) may be computed in the high tem- 
perature phase by a calculation parallel to that of [ID]. On the constrained 
ring a perturbation n a (x) = 1/3 + a a cos(2irmx) + b a sin(27rmx) of the con- 
stant solution satisfies (|3.4p and ^ n a (x) = 1 if and only if m = 3k + j for 
j = 1 or 2, and 

a a+ i = --a a + —b a , b a+1 = --b a - (-l) j —a a . (4.12) 

One may thus treat cla and 6^ as the independent parameters. The proba- 
bility of the profile {n a (x)} is exp{— SNJ 7 ^ 3 ^ ({n})}, and to quadratic order 
in the perturbation, 

T m ({n\) ~ constant + \27rm + (-lV'V3/3l (t& + (4.13) 

87rm L 



Thus 



la\) = (b\) = — — , (a A b A ) = 0. (4.14) 

27iY(2vrm + (-1^3/3) X 7 V 7 



Summing over all the fluctuations, i.e., over m, we obtain 
(n a (x)n a (y)) c = f y ^ cos[2nm(x - y)] 



(4.15) 



m=3k+j 



All connected two-point functions (n a (x)n 1 (y)) c may be obtained on the 
constrained ring from ()4.15j) via (j3.4j) , and then on the interval using (|4.1ip . 
Note that KTEb diverges as (3 / j3 c . 

4.1 The canonical free energy F(r) 

The free energy in the canonical model, for mean densities r A , rs, rc satis- 
fying ta + tb + tc = 1, is given by 

F(r) = F(r A ,rB,r c ) = min T({n(x)}), (4-16) 

{n(x)} 
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with the minimum taken over all profiles n{x) satisfying the constraints 
(|2.5p . The grand canonical free energy may then be computed in two ways: 

F(X) = inf ' f({n}) (4.17) 

{n(x)} 




F(r)-J2^r a \, (4.18) 



where the infimum in (|4.17p is over all profiles. We can obtain information 
on the structure of F(r) from the above results for the minimization problem 
(|4.17p . together with the trivial remarks that a unique minimum for (|4.17p 
implies a unique minimum for (14.18P and that such a unique minimum im- 
plies that the surface y = F(r) lies above the plane y = F(X) + ^ a X a r a 
and touches it at a single point. 

In particular, the fact that when f3 < (3 C there is for all A a unique 
minimizer for (I4.17P implies that for such (3 the function F(r) is convex. 
When f3 > f3 c the minimizer for (|4.17j) is unique except in the case Xa = 
Xb = Ac* = 0, when the plane mentioned above is horizontal. In that case 
the minimum occurs at points lying above a certain simple closed curve 
T (= Tp) in the plane Yl a Ta = ^' w hh the point ta = tb = tq = 1/3 in 
its interior; sample curves are shown in Figure EJ V may be parametrized 
as r*(t), < t < tk, where K is the parameter in the type 1 solution of 
Theorem |4JJa) and 



o rr K /6 

r* a (t) = — y K (s + t a + t)ds. (4.19) 

TK J-t k /6 



(The fact that this curve is simple follows, for example, from Propo- 
sition lA.lf dl.) The three- fold symmetry then implies that the surface 
y = F(r) has a "tricorn" shape. 



5 Concluding remarks 

It is natural to compare the phase diagram obtained here for the one di- 
mensional reflection asymmetric ABC model with that of the correspond- 
ing symmetric model, that is, the mean field three state Potts model (see 
[T].|l7|). We will define the latter by replacing the sum over j > i in (jl.ip 
by a sum over all j / i and dividing by 2; this yields 

^(0 = t^WaNc+NcNb+NbNa] = -±-[N 2 -{N 2 A +Nl + N 2 c )}. (5.1) 
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'A 



Figure 2: Curves T$ in the r^-rs plane along which F(r) achieves its mini- 
mum value, for ft = 3.75, ft = 4-25, and ft = 6.05 (ft = 2ir/^3 ~ 3.63). 



The energy (|5. 1|) is related to that of the standard mean field Potts model 
[17] by a choice of energy scale and a shift of the ground state energy. It is, 
as is usual for mean field models, independent of dimension and geometry. 
There is thus no spatial structure in the system and the canonical measure 
just gives equal weights to all configurations. 

The canonical free energy functional with prescribed values of r a — 
Jq 1 n a (x) dx is 

F*({ n }) = ^[tatc + r B TA + revs] - 5({n}), (5.2) 

with S ({n}) still given by (|2.4p . For all (3 the minimizers of J-* are the 
constant density profiles p a (x) = r a , and there are no phase transitions of 
any kind in the canonical system. The corresponding minimum value 

a a a ct 

of J-"*({ra}) is in fact just the value of J-"({n}) evaluated at these constant 
profiles (this follows from our choice of the factor 1/2 in (|5.ip ). and hence 
is an upper bound for the free energy F(r) of (|4.16p . 

The situation is quite different for the grand canonical ensemble. Here 
the analogue of (|4.17p is 

F*(A) = inf{F*(r)-^A a r a } (5.4) 

a 
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The analysis of F*(X) leads to a phase diagram completely different from 
that of the reflection asymmetric grand canonical model considered in Sec- 
tions 3 and 4 above pQ. In particular, (|5.4p exhibits a first order phase 
transition for A^ = Xb = Ac at f3* = 8 log 2. For /3 < (3* the minimizer is 
r A = r B = r c = 1/3; for {3 > f3* there are three minimizers, each rich in one 
of the three species, and at f3 = (3* all four of these states are minimizers. 

5.1 Higher dimensions 

As was already noted and is well known, the standard mean field models 
with symmetric interactions do not depend on the dimension or topology of 
the spatial structure of the system considered. This is clearly not the case for 
models with reflection asymmetric interactions, such as the one-dimensional 
ABC model considered in this paper. 

We comment now on various possible generalizations of such reflection 
asymmetric mean field models to higher dimensions, taking for simplicity 
the dimension to be two and the lattice to be an iV x JV square in 1? . Let us 
consider first a situation in which the mean field interactions are symmetric 
in the vertical direction but of the form (jl.ip in the horizontal direction. 
This yields an energy of the form 

N N-l N 

%) = pEEEE^^ao (5.5) 

k,l i=l j=i+l ol 
N-l N 

= Yl Yt^^Va+lij), (5.6) 

i=l j=i+l a 

where 

1 N 

Va(i) = ^2va(i,k). (5.7) 

k=l 

The energy functional £ obtained from (|5,6p in the scaling limit is identical 
to that given in f)2.3f) with n a {x) replaced by n a {x) = J Q n(x,y) dy. The 
entropy term (compare (|2.4p ). 

-S = y2 / n a (x,y)logn a (x,y)dxdy, (5.8) 

is clearly minimized, subject to a specified {n a (x)}, by setting n a (x,y) = 
n a (x), and so density profiles which minimize f3£ — S depend only on x and 
are the same as for the one dimensional case, both for the canonical and 
grand canonical ensembles. 
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Remark 5.1. The two-dimensional chiral clock model [21 S] also contains 
interactions — in that case, nearest-neighbor ones — which are reflection sym- 
metric in the vertical direction but not in the horizontal one. When the 
parameter A (in the notation of [2]) has value 1/2 the energy, up to an 
additive constant and a rescaling, is 

N-l N N AT-1 

^2 ^2^2va(i,k)r] a+ 2(i + l,k) - ^2^2 ^2va(i,k)r] a (i, k + 1), (5.9) 

i=l k=l a i=l k=l a 

so that the interactions in the horizontal direction have a form reminiscent 
of ([LTD. 

A second possibility is to take the reflection asymmetry to be the same 
in the x and y directions. In this case (jl.ip takes the form 

^ N-l N 

E N 2 (C) = ^EE Yl r, a (i,k)r, a+2 (j,l). (5.10) 

a ik=l J=*+l 
l=k + l 

The analysis of this model seems considerably more complicated and we will 
attempt no discussion here. 
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A Proof of Theorem 14.1( b) 

We begin by giving an alternate form of the boundary conditions (14. 7|) . With 
(p^|) and (j2~T0j) these become 



K - K+i = log{y K (t a - /3)) - log(y K (t a + f3 - t k /3)) . (A.l) 
From ([221) and flZgJ, (logy K (t))' = [y K (t + T K /3)-y K (t-T K /3)]/2, so that 

K ~ K+i = I [ \vk (t + ^f) - y K (t - ?f)] dt. (A.2) 



2 



t a +P~T K /3 



3 / V 3 



The solution of (|A.2[) which also satisfies X a = is 

I fs a + {T K /6-l3) /l \ 

A a = - / o - Vi<(t) dt, (A.3) 



2 ■ls a -(T K /G-P) V 3 
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where s a = t Q + tk/2. The form ()A.3|) is convenient when 2/3 < tk/3; if 
2/3 > t/^/3 we may rewrite this as 

A Q = -/ [VK(t)--)dt. (A.4) 



2 7 Sq _ (/3 _ ta ,/ 6 ) V 3 , 

The representations (lA.3j) and (lA.4j) are useful because they translate the 
boundary conditions for the grand canonical model into a form similar to 
the condition (|2.1ip in the canonical model. 

We need also to recall from p3] some further properties of the function 
Uxit) and its definite integrals 

Y(K, s, 5) = f S+S y K (t) dt and W(K, s, S) = f ^ y K (t + ^) dt. 

J s- — 5 J s — S 

(A.5) 

Note that from (pTTT]) . 
and that from (f2TT0|) . 



1 Y(K,t a ,j3) (A.6) 



2/3 



K = ^-^Y(K,s a ,5) = ^-^W(K,s a+l ,5), (A.7) 
for 5 = t k /6 - > 0, while for 6' = -5 > 0, 

A* = ^Y(K,s a ,5')- 6 - = ^W(K,s a+1 ,S')- d -. (A.8) 
Proposition A.l. For < K < 1/27: 

(a) (i) yxit) is even and tk -periodic (and hence also symmetric about t = 
tk/2), takes its minimum value at t = 0, is strictly increasing on [0,tk/2], 
and takes its maximum value at t = tk/2. Moreover, (ii) yx(t — t k/3) + 
VK(t) + VK(t + Tfc/3) = 1 for all t. 

(b) The minimum value a = a(K) = yx(0) of yx is an increasing function 
of K satisfying < a(K) < 1/3. The maximum value b = b(K) = yK^jKl^) 
isb= [2-a-V4a-3a 2 ]/2 ; andy K (T/6) = (l-6)/2, Vk(t/3) = (l-a)/2. 

(c) (i) For fixed K and 5, with < 5 < tk/2, the function Y(K,t,8) shares 
with VKif) the properties listed in (a.i). Moreover, (ii) 

Y(K, t - t k /3, 5) + Y(K, t, 5) + Y(K, t + r fc /3, 5) = 26. (A.9) 

(d) For < 5 < Tk/2, Y(K, 0, 5) is strictly decreasing, and W(K,Tk/6,8) 
strictly increasing, in K. 
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Finally, for < K 2 < K x < 1/27: 

(e) (i) For any to the curves yidito + t) and yK 2 (t) intersect exactly once in 
the interval < t < tk 2 /2, and (ii) yK 2 (t) < UKi(t) and yK 2 {jK 2 /2 ~ t) > 
y Kl (T Kl /2-t) forO<t<T Kl /6. 

Proof. These results either appear in p3] or are immediate consequences 
of results appearing there. For (a) and (b) see Section 5.2 of |14| and in 
particular Remark 5.1(a); for (c.i) see Remark 5.3(b). (c.ii) follows from 
(a.h). The first statement of (d) follows from the fact that Y(K,0,5) is 
continuous in K and, for < 5 < Tfc/2, approaches 25/3 as K f~ 1/27 
and as K \ 0, together with Theorem 6.1 of p3] which, if one takes 
there ta = re, asserts that for given r# with < r# < 1/3 there is 
at most one value of K satisfying (|A.6p . The second statement of (d) is 
verified similarly. Finally, (e.i) is a special case of Lemma 6.2(a) of |14j 
and (e.ii) then follows from (e.i) and the inequalities yK 2 (0) < 2/^(0), 
Vk 2 {tkJQ) < yK 2 (TK 2 /G) < yKi(r Kl /Q), yK 2 {TR 2 /2) > y Kl { T K 1 /2), and 
yK 2 {TR 2 /2 - r Kl /6) > y K2 {TK 2 /3) > 2/^(^/3), easily obtained from the 
properties given in (a) and (b). □ 

We now turn to the proof of Theorem 14.1( b). We know (see the remarks 
at the beginning of Section H|) that there is at least one minimizer and that 
every minimizer satisfies the ELE (|4.6p . (|4.7p . Thus the conclusion of the 
theorem will follow from: 

Lemma A. 2. If \a, ^b, and Ac are not all equal then: 

(a) No solution of (|4.6[) . (14. 7|) of type n, n > 2, can minimize T . 

(h) At most one solution of (|4.6p . (j4.7j) of type 1 exists. 

Remark A. 3. In proving Lemma IA.2I we need not consider either the 
constant solution of the ELE or nonconstant solutions for which 2/3 = 
(n — 2)/3 c /3, both of which satisfy (|4.7|) only when all the X a are equal. 
We may also suppose, without loss of generality, that 

Aj4 < Ac < Xb, with Xa < Ac or Ac < Xb- (A. 10) 

If S = t k /6 - f3 > it then follows from ([A/Tl) that Y(K,s B ,S) < 
Y(K,sc,S) < Y(K,sa,S) and then from Proposition lA.l( c.i) (see FigureEl 
which displays graphically the qualitative properties of Y(K,s,5) implied 
there) that < sb < t~k/6, so that tk/2 < ts < 2tr-/3 and hence, from 
Proposition lA.lf c.i) and (IA.6p . that < rc < tb- If A^ < Ac then 
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Figure 3: Plots showing qualitative features of Y(K,s,S) (solid), Y(K,s + 
tk/3, 8) (dotted), and Y(K,s — tk/3,8) (dashed) for < 5 < tk/2, based 
on Proposition I A , 1 f c . 1 ) . 

sb > 0, ts > t~k/2, and ta < rc; similarly rc < re, if Ac < \b- Simi- 
larly, if tk/6 < /3 < tk/2 then (now using (jA.8p ) tk/2 < sb < 2tk/3 and 
rB<rc<rA-, again with strict inequality for two of the X a implying the 
corresponding inequality for the r a . 

Proof of Lemma I A. 2( a): Consider some type n solution p(x), n > 2, of 
fljgP , (l47j) : p(x) has the form t$2l$) with 2/3 > t k /3. We need to find a 
profile p(a?) with .F({p}) < .F({p}). There are three subcases: 
Case (a.i) 2/3 > tk- In this case it was shown in p3] that there is a 
rearrangement p(x) of p(x) with ^({p}) < .F({p})- This rearrangement 
does not change the mean densities r a and hence also F{{p}) < T({p}). 

Case (a.ii) 2/3 = tk- In this case the solution p(x) has mean densities 
r a = 1/3, so that ^A Q r a = 0. From the description of the curve V in 
Section 14.11 it follows that for some z > there exists a minimizer p(x) of 
P^^ with mean densities f a = 1/3 + zA Q , so that ^ A Q f a > 0. But then 

^ X \{ P }) = tw> H{p}) > f^°H{p}) > f^' X} ({p}). (a.ii) 

Case (a.iii) tk > 2/3 > tk/3. By Remark IA.31 r\e < rc < r^, with 
< if Ac* < As and rc < if A a < Ac- Consider now the profile p with 
Pa (#) = Pa+i( x )- The canonical free energy functional satisfies F{{p}) = 
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and so 

Hip}) ~ £({p}) = Y Xara ~ Y A « r °+i 

a a 

= (Xa ~ Xc)(r A ~ r B ) + (A B - \ C )(r B - r c ) < 0. (A.12) 

□ 

The next result, the key to the proof of Lemma IA.2f b). gives certain 
monotonicity properties of Y(K, s, 5) and W(K, s, 5). 

Lemma A. 4. If K, s, and 5 satisfy < K < 1/27 and < s,5 < tk/6, 
then: 

(a) For fixed K and s the function 25/3 — Y(K,s,5) (respectively 25/3 — 
W(K,s,5)) is strictly increasing (respectively strictly decreasing) in 5; 

(b) For fixed K and 5 the functions Y(K,s,5) and W(K, s,5) are strictly 
increasing in s; 

(c) For fixed s and 5 the function Y(K,s,5) (respectively W(K, s, 5)) is 
strictly increasing (respectively strictly decreasing) in K. 

Proof, (a) We rely throughout on Proposition lA.l( a.b), From < s + 5 < 
tk/3 and — < s — 5 < tk/G it follows that vk{s + 5) < (1 — a)/2 and 
y K (s -<&)<(!- b)/2. Then from ijAlH) . 



d_ 

d5 



^-Y(K,s,5) 



\-yK{s + 5)-y K {s-5)>^- l ->Q, (A.13) 



as is easily verified from b = [2 — a — \/4a — 3a 2 ]/2 with < a < 1/3. To 
show that (d/d5)(25/3 — W(K,s,5)) < it suffices similarly to verify that 

z(K, s, 5):=y K ( S + ^- + 5)+y K (.s + ^--5)> | (A.14) 

Because yx is even and r^-periodic, z is invariant under (s,5) — > (s',5') 
with s' = tk/6 — 5, 5' = tk/6 — s, so that it suffices to verify (|A.14p for 
s + 5 < tk/6, and since under this condition both terms in z(K,s,5) are 
increasing in s it suffices to consider s = 0. But because yx is even, 



z(K,0,5) -- o 



Vk[ — + + VK y — — + 5 

+VK (--ir -o) +VK ( ~!f - 5 



3 J " \ 3 
1 " \ [VK(S) + y K (S)] > ^ > 2/3- (A.15) 
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(b) See Proposition lA.lf c). 

(c) The proofs for Y and of W are similar and we check only Y. Suppose 
that < K 2 < K\ < 1/27 and that for some s* E [0,77^/6], 



Y(K!,s*,5) < Y{K 2 ,s*,8). 



(A.16) 



Then certainly yKi(t*) < J/K 2 (**) f° r some t* £ [s* — 6, s* + 5], and since 
y.?fi(0) > yx 2 (0)' Proposition lA.lf e.i) implies that UK^t) < VK 2 {t) f° r * 6 
Then for s G [s*, 77^/6], 



as 

= [VKi (s-S)- y K2 (s - 5)] + [vk 2 (s + 8)- DKi (s + 5)} 

> 0, (A.17) 

since both terms on the right had side are positive. But for < 5 < 77^/6, 

W{K x ,t k JS, S) < W(K 2 ,tkJ6,6), (A.18) 

by Proposition IA.lf e.ii). and so from Proposition lA.lf a). 



Y(Ki,t Ki /6,S) 



1 r T K 1 /Q+5 
1 /" r A' 1 /6+5 r 



VKi ( * - ) + VKi (t) 



dt 



1 - .V/v , [ t+ T -f 



dt 



TK X /6— S 

> 5-Uv(K 2 ,t k JQ,5) 
= Y(K 2 ,T K2 /6,6) 

> y(^ 2 ,r^/6,(5) (A. 19) 
since 77^ < tt< 2 (see Remark l2.ip . contradicting ()A. 16|) and (|A.17p . □ 

Proof of Lemma \A. 13( b): For type 1 solutions we have from (|A.7p that 

5 1 



Xa = ---W(K,s b ,S), 



(A.20) 



with 5 = 77-/6 — (3 > and, by Remark IA.31 < sb < TR-/6. Thus 
the existence for some A of two type 1 solutions would correspond to the 
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existence of (Jfi,si) and (K 2 ,s 2 ) with < K 2 < K\ < 1/27 and < < 
TK J6, i = 1, 2, such that 2<5i/3 - K(Jfi, si, £1) = 2J 2 /3 - Y(K 2 , s 2 , 5 2 ) and 
2<5i/3 - W(ETi,si,5i) = 25 2 /3 - W(if 2 , s 2 , <5 2 ), where S { = tjq/6 - /3 for 
2 = 1,2. Then from Lemma I A.4( a.c) . 

^-Y(^ 2 , S2 ,5 2 ) = Ml-y^,^) 

< ^-Y(K 2 , S1 ,5 2 ), (A.21) 
so that Lemma lA.4( b) implies that s% < s 2 . But also 

^-W(K 2 ,s 2 ,5 2 ) = ^-W(K 1 ,s 1 ,5 1 ) 

> ^-W(K 2 ,s 1 ,S 1 ) 

> ^- W (K 2 , S1 ,5 2 ), (A.22) 
implying that s\ > s 2 , a contradiction. □ 
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